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Method and Apparatus for Improving Transmit Antenna Weight Tracking Using Channel 
Correlations in a Wireless Communication System 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application is a continuation-in-part of U.S. Application No. 09/632,081, filed 
08/02/2000, entitled "Adaptive Antenna Method and Apparatus", hereby incorporated 
by reference herein for its teachings on digital communication systems. This 
application claims the benefit of the priority filing date of U.S. Application No.: 
09/632,081 in accordance with 35 USC § 120. This application is related to and also 
claims the benefit of U.S. Provisional Application No. 60/278,501, filed 03/23/01, 
entitled "Adaptive Antenna System and Method". 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

This invention relates to digital wireless communication systems, and more 
particularly to methods for improving transmit antenna weight tracking using channel 
correlations in a wireless communication system. 

2. Description of Related Art 

Digital wireless systems presently being deployed are commonly seen to be forward 
link (i.e., base station ("BS")-to-mobile station ("MS")) capacity limited. This is 
mostly due to asymmetric traffic scenarios: most data applications, such as web 
browsing, are envisioned to be sending relatively large quantities of data to the mobile 
user on the forward link, while the reverse link may contain only small control packets 
(i.e., a new IP address to download) so that the forward data rate is much greater than 
the reverse data rate. 
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Forward capacity can be increased with an array of antennae in a base station whose 
individual carrier amplitudes and phases can be adjusted based on feedback from a 
MS. In particular, a transmit adaptive antenna (TxAA) algorithm will typically 
transmit the same MS-specific waveform on multiple antennae, applying adaptive 
complex weights to the baseband signal applied to each antenna. To allow the MS to 
perform coherent demodulation, a dedicated pilot channel is typically transmitted in 
the same manner as the data. 

While receive adaptive antenna (RxAA) algorithms at the BS are fairly straight 
forward, TxAA algorithms are not. The mere definition of an "optimal" TxAA 
algorithm is not unambiguous as the optimization of one mobile's forward link can 
degrade another's, leading to complex tradeoffs which are not part of the RxAA 
problem. 

Most TxAA algorithms will require some a priori knowledge of the transmission 
channel from the BS to the MS. Given this knowledge and a defined "optimality" 
criterion, the BS can determine the TxAA weights. The difficulty is that, although the 
MS can measure the channel, it is the BS which needs the measured information to 
adjust its transmit weights. Also, there can be some additional complexity in that the 
MS must generally measure the channel of each transmit antenna separately, in 
addition to measuring the channel of the overall transmit (Tx) weight-adjusted signal. 
The latter measurement is required for the demodulation of the signal. 

A number of methods for implementing TxAA algorithms have been proposed. 

According to some of these proposals, a few bits are allocated to the MS to encode the 

channel gain and phase. See, for example, Thomas Derryberry, Balaji Raghothaman 

(Nokia) "Transmit Adaptive Arrays without User Specific Pilot", document # C30- 

19990817-030, submitted to 3GPP2 August 1999; Mark Harrison (Motorola) "Tx AA 

Parameter Recommendations", document #C30-19990914-010, submitted to 3GPP2, 

Tokyo Japan, September 1999; and Mark Harrison, Kiran Kuchi (Motorola) "Open 

and closed loop transmit diversity at high data rates on 2 and 4 elements", document 

#C30-19990817-17, submitted to 3GPP2 August 1999. These methods do not allow 
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the desired antenna weights to be precisely determined because the channel state must 
be distorted in order to be fed back to the BS with a low bit rate. 

Typically, the academic literature has assumed that the full channel information is 
available at the transmitter, which is not a practical assumption. See, e.g., Jen- Wei 
Liang, Arogyaswami Paulraj "Forward link antenna diversity using feedback for 
indoor communication systems" Proceedings, 1995 International Conference on 
Acoustics, Speech, and Signal Processing, May 1995; Farrokh Rashid-Farrokhi, K.J. 
Ray Liu, Leandros Tassiulas "Transmit beamforming and power control for cellular 
wireless systems" IEEE Journal on Selected Areas in Communications, Vol. 16, No. 
8, October 1998. There have been several submissions to the TIA standardization 
body for TxAA algorithms with channel feedback from MS to BS, some of which are 
referenced above. All of these submissions require that the MS measures a primary 
and secondary pilot transmitted from the primary and secondary antennae. Some 
submissions have shown that 4 antenna transmission can yield improved performance, 
however, with no recognition of the increased complexity at the MS. 

The use of a few bits of channel feedback leads to substantial degradation relative to 
possible performance, because the feedback requires low bit rate quantization. Also, 
the MS must individually measure and transmit information for each antenna; the MS 
requires extra hardware to perform these functions. Furthermore, the above-described 
methods do not gracefully grow to accommodate an increase in the number of 
antennae. 

The above-mentioned systems are illustrated in FIGURES 1-3, which illustrate a 

CDMA system. FIGURE 1 shows a transmitter 10 with two antennae, a first antenna 

(antenna #0) 12 and a second antenna 14 (antenna #1). As is shown, different 

common pilot signals are associated with different antennae. For each user, such as 

user #0 and user #1, the transmitter 10 includes an adder, e.g., the adders, 16 and 18 

respectively, that adds together a dedicated pilot signal for the user and the forward 

traffic for that user. For user #0, multipliers 20 and 22 multiply the summed signal 

with complex weights for the first and second antennae 12 and 14, respectively. For 
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each antenna, the weighted user signals are summed by adders 24 and 26 and the 
result is added by adders 28 and 30 to the pilot signal for that antenna. It should be 
noted that FIGURE 1 represents a complex baseband equivalent, as no RF modulation 
stage is shown. 

5 

FIGURES 2 and 3 show possible embodiments for receivers that may be used to 
receive signals transmitted by the transmitter 10. As shown in FIGURE 2, a received 
signal is divided into three components: one signal corresponding to the dedicated 
channel, one signal corresponding to the common pilot for the first antenna and one 

10 signal corresponding to the common pilot for the second antenna. This division is 
accomplished by multipliers 32, 34 and 36 and accumulators 38, 40 and 42. A more 
hardware efficient embodiment is shown in FIGURE 3, which employs multiplexer 44 
and demultiplexer 46 to alternately select between the signals for the different transmit 
antennae. The time-multiplexed processing used by the embodiment of FIGURE 3 

15 saves hardware at the expense of 3dB loss of precision for each channel. This tradeoff 
may be acceptable when a low bit rate, low precision channel reporting is used by the 
mobile to report these channel estimates to the BS. 

The channel estimate attained by the mobile is coded into a low bit representation. 
20 The bit rates mentioned in the above cited references are 1, 2 or 4 bits (lb phase, 2b 
phase, or 3b phase + lb amplitude). So, for a forward channel vector c, the mobile 
generates the estimate c , which is then quantized to produce the feedback estimate I . 

Note that the MS channel estimation hardware of FIGURES 2 and 3 can be used for 
25 an arbitrarily large number of Tx antennae, provided that there are unique pilot codes 
for each antenna and that the MS has information regarding these codes. 

Finally, the proposed systems have the MS report the channel estimate based on 1 

"path". In the presence of resolvable multi-path due to delayed reflections of the 

3 0 transmitted waveform, particularly for CDMA, there may be more than one path that 

is usable to the MS. In order to report the channel for N such paths, the MS-to-BS 

feedback rate would be required to increase N-fold, and the number of such paths 
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would somehow need to be communicated to the BS. This solution is not practical 
however, and instead, the MS reports the channel estimate for only the strongest path. 
This solution discards some useful channel characteristics which could further 
increase performance under these circumstances. 

5 

The algorithm employed by the BS to utilize the received channel information would 
most likely be a simple matched transmission weighting. That is, the forward weights 
chosen would be the conjugate of the forward vector channel, so that the weights are 
determined as follows: 

10 w = c. 

This formulation attempts to maximize signal power to the mobile without regard to 
□ the locations of other mobiles. It does not steer nulls to the other mobiles. The 

jig channel estimate from a given mobile could be used to determine transmission nulls of 

!;3 15 other mobiles, but the coarseness of the channel estimate (no greater than 4 bits) 

I'll makes this approach ineffective. 

CS 

!;H Schemes such as that shown in FIGURES 1-3 have significant disadvantages. In 

III 

fU particular, as previously mentioned, the use of a few bits of channel feedback requires 

j^l 20 quantization, which leads to substantial degradation relative to the possible 
HO performance. Also, because an MS unit must individually measure and transmit 

information for each antenna, the MS requires extra hardware to perform these 
functions. Furthermore, the above-described methods do not grow gracefully in 
application to more antennae. 

25 

One method used by the present invention in implementing a TxAA algorithm 
(referred to as the "Adaptive Antenna Method and Apparatus 11 ) that overcomes these 
disadvantages is described in the above-incorporated U.S. Provisional Application No. 
60/278,501, U.S. Application No. 09/632,081 and hereinbelow. In accordance with 
3 0 this method, a transmitter includes hardware and/or software for adaptively updating 

weights for a plurality of antennae. In particular, according to the present invention, 
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the transmitter generates a pair of test weight vectors, each vector comprising a 
plurality of complex entries, with each entry corresponding to a different one of a 

yy _|_ ft.y 

plurality of antennae. The first complex weight vector is equal to ,, base » and 

K-+>M 

yy — R , y 

the second complex weight vector is equal to .. base — f , where v is a test 

Whose- PA 

perturbation vector comprising a plurality of complex entries, with each entry 
corresponding to a different one of a plurality of antennae, p is an algorithm constant 
and w base is a vector that is updated based on feedback received from a receiver. The 
first and second complex weight vectors are applied to a dedicated pilot signal during 
alternate time intervals. During each time interval, the average of the first and second 
complex weight vectors is applied to the data traffic transmitted by the transmitter. 

A receiver (e.g., a mobile station or ("MS")) alternately receives the pilot signal as 
multiplied by the first and second weight vectors as described above. The mobile 
station determines which of the weighted pilot signals resulted in a stronger signal 
received at the mobile and, based upon this determination, transmits feedback 
information. The transmitter receives the feedback and updates the first and second 
weights accordingly. In particular, if the first weight resulted in a stronger signal, w base 

w + B • v 

is updated to become „ base — - — Tx (the previous first weight) whereas if the second 

yy — R . y 

weight resulted in a stronger signal, w is updated to become „ base j, (the 

previous second weight). The above-described process is repeated with the new w base . 
New test vectors v are generated and applied after each MS channel measurement is 
made. 

A disadvantage of the adaptive antenna method and apparatus and the other TxAA 

techniques described above with reference to FIGURES 1-3 is that these techniques 

may track weight vectors at a relatively slow rate. These techniques provide no 

mechanism for using partial knowledge of the channel, such as the channel gain vector 

autocorrelation. Thus, a need exists for a method and apparatus for improving transmit 
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antenna weight tracking. Such a method should utilize existing TxAA techniques to 
improve the performance of a wireless communication system. The inventive method 
should not degrade performance due to quantization from utilizing a few feedback 
bits. In addition, the method and apparatus should be able to accommodate multiple 
antennas. 

The present invention provides such a method and apparatus for improving transmit 
antenna weight tracking. 
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SUMMARY OF THE INVENTION 

The present channel vector autocorrelation tracking technique improves weight vector 
tracking in a wireless communication system by taking advantage of knowledge of the 
correlation between elements of the channel gain vector. Several embodiments of the 
present invention are described hereinbelow. 

In one embodiment of the autocorrelation tracking technique, the present invention 
extracts a biased gradient estimate by utilizing a correlated perturbation vector 
probing technique, where the perturbation vector is generated with a autocorrelation 
derived from an estimate of the channel vector autocorrelation matrix. The 
autocorrelation with which the perturbation vector is generated is selected such that 
the stronger subspaces of the channel autocorrelation are emphasized in the 
perturbation probing and update. This is desirable because these subspaces are more 
likely to be occupied by the channel vector at any given time than weaker subspaces, 
and thus the gradient estimate is biased to these emphasized subspaces. Specific 
embodiments use a perturbation autocorrelation given by an exponentiation of the 
estimate of the channel matrix added to an equi-diagonal matrix (i.e., scaled identity 
matrix). 

The details of the embodiments of the present invention are set forth in the 
accompanying drawings and the description below. Once the details of the invention 
are known, numerous additional innovations and changes will become obvious to one 
skilled in the art. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE 1 is a simplified block diagram of one type of conventional two-antenna 
transmitter that employs antenna weighting derived from receiver feedback 

FIGURE 2 is a simplified block diagram of a first type of conventional receiver that 
estimates the channel for each transmit antenna individually. 

FIGURE 3 is a simplified block diagram of a second type of conventional receiver 
that estimates the channel for each transmit antenna individually. 

FIGURE 4 is a simplified block diagram of one embodiment of a transmitter 
constructed according to the present invention. 

FIGURE 5 is a simplified block diagram of one embodiment of a receiver constructed 
according to the present invention. 

FIGURE 6 is a block diagram of a preferred embodiment of a receiver constructed 
according to the present invention. 

FIGURE 7 is a flow chart that outlines pertinent features of the operation of the 
transmitter shown in FIGURE 4. 

FIGURE 8 is a flow chart that outlines pertinent features of the operation of the 
receiver shown in FIGURE 5 

FIGURE 9 is a flow chart that outlines pertinent features of the operation of the 
receiver shown in FIGURE 6. 

FIGURE 10 is a timing diagram that shows a preferred timeline for weight 
adjustment, channel measurement and feedback. 
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FIGURE 1 1 is a block diagram of an exemplary device that is adapted for use with the 
present invention. 



Like reference numbers and designations in the various drawings indicate like 
elements. 
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DETAILED DESCRIPTION OF THE INVENTION 

Throughout this description, the preferred embodiment and examples shown should be 
considered as exemplars, rather than as limitations on the present invention. 

The method and apparatus of the present invention improves transmission weight 
vector tracking in a wireless communication system. The present invention improves 
weight vector tracking by utilizing one of several techniques described hereinbelow. 
The present invention can utilize one technique or any combination of techniques to 
improve weight vector tracking. According to this invention, an estimate of the 
autocorrelation matrix of the transmission channel gain vector is generated, and the 
perturbation vectors used for probing and gradient feedback are generated with an 
autocorrelation which is a function of the channel autocorrelation matrix estimate. 

The present inventive method and apparatus is primarily intended for use in digital 
wireless communication systems. An exemplary DS-CDMA communication system 
adapted for use in the present generating transmit adaptive antenna weights using 
binary gradient feedback method and apparatus is described. However, as those of 
ordinary skill in the art shall recognize, any convenient wireless communication 
system can be used to practice the present invention. The communication system 
should allow the mobile station the capability of measuring radio frequency signal 
strengths and providing binary feedback to the base station. The present invention can 
be used in any wireless communication system where improving the tracking of 
weight vectors is desirable. 

The following table defines some terms used throughout this description of the 
invention. 
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Nomenclature: 

Vector and matrix variables are denoted in bold, scalars are not in bold. 
Superscript "T" represents matrix or vector transposition. 

Superscript "H" represents matrix or vector hermitian transpose (or conjugate 



transpose) 



Term 


Definition of Term 


AA 


Adaptive antenna, more generally adaptive antenna aigoninm. 


BS 


Base Station, fixed station which communicates with multiple 
mobile stations and provides their gateway to the fixed network. 


CDMA 


Code Division Multiple Access. 


DS-CDMA 


Direct Sequence CDMA. 


Forward link 


Radio link from the transmitting base station to the receiving 
mobile station. 


HW 


Hardware, referring to the fixed digital logic portion of a modem 
(as opposed to SW). 


MS 


Mobile Station, a user of the cellular network, e.g. a cell phone or a 
wireless data terminal. 


pilot 


A known signal or "training sequence" transmitted so that the 
receiver may use it to estimate the channel response. 


reverse link 


Radio link from the transmitting mobile station to the receiving 
base station. 


Rx 


Receive. 


SW 


Software, typically to distinguish the (modifiable) SW portion of a 
modem from the HW. 


Tx 


Transmit. 



Transmit Adaptive Antenna Weightin g Embodiment 

For clarity purposes, certain non-essential assumptions are made in the present 
detailed description. First, the wireless communication system described is assumed 
to be a DS-CDMA (Direct Sequence, Code Division Multiple Access) system that 
uses "pilots" (or "pilot channels") for use in performing coherent demodulation. In 
fact, most modern digital wireless systems include a pilot channel (sometimes known 
as a "training sequence", when the pilot channel is time multiplexed, rather than code 
multiplexed) which is used for channel estimation. Any such system can be used in 
practicing the present invention. 
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Second, the below described adaptive array transmitter is considered to exist in a base 
station and the receiver is a mobile station. This applies to cellular communication 
applications. In one embodiment, the mobile station may also have a transmit array 
operating in the same manner. Alternatively, the system that uses the present 
inventive method and apparatus may not be a cellular system, but rather, it may be a 
peer-to-peer communication system such as those defined by the Bluetooth™ 
Specification. 

The invention as described can be applied without regard to how many antennae the 
transmitting BS is using, and in particular the MS does not need knowledge of the 
number of antennae being used. 

FIGURE 4 is a simplified block diagram of a BS transmitter 50 mode according to one 
possible embodiment of the present invention. For ease of illustration, two antennae, 
a first antenna 52 and a second antenna 54 are shown, but it will be appreciated that 
the present invention may easily be scaled to any number of antennae. Also for ease 
of illustration, although BS hardware corresponding to one mobile unit is shown, the 
transmitter 50 will typically have analogous hardware for a number of other users. 

As shown in FIGURE 4, the transmitter 50 comprises a plurality of multipliers 56, 58, 
60, 62, 64 and 66, summers 72 and 74, and multiplexers 68 and 70. Multipliers 64 
and 66 multiply a dedicated pilot signal and a forward traffic signal, respectively, with 
corresponding gains ("real valued"). Multipliers 60 and 62 multiply the forward 
traffic signal by weights that are components of the weight vector w, w 0 and w„ where 
w 0 is the weight for the first antenna 52 and Wi is the weight for the second antenna 54. 
The weight values are described below in more detail. w<, !0dd and w 0<even , weights from 
the weight vectors w odd and w even respectively, are provided as inputs to the 
multiplexer 68, which alternates between them. Similarly, w Uodd and w 1)CTen weights are 
provided as inputs to the multiplexer 70, which alternates between them. The 
weighted traffic and pilot signals are summed by the summers 72 and 74 and 
transmitted by the first antenna 52 and the second antenna 54, respectively. The 
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mathematical representation of the waveform transmitted by the antennae 52 and 54 is 
described below in more detail. 

FIGURE 5 is a simplified block diagram of an MS receiver 80 according to one 
possible embodiment of the present invention which demonstrates the principle of the 
feedback decision made by the MS. The receiver 80 comprises an antenna 82, a 
multiplier 84, and an accumulator 86. The receiver 80 further comprises a 
demultiplexer 88, accumulators 90 and 92 and a summer 94, which are preferably 
implemented in software. Signals received by the antenna 82 are multiplied by a pilot 
despreading code, as is well known in DS-CDMA systems, and the result is passed to 
accumulator 86, which accumulates a symbol from the chips it receives. "Even" and 
"odd" channel estimates are alternately accumulated in the accumulators 90 and 92. 
The even and odd estimates (i.e., the estimates for the alternating even and odd time 
slots) are summed by the summer 94 to generate a dedicated channel estimate (for 
coherent demodulation). Also, in a decision block 96 the power of the even and odd 
estimates are compared and the value of a feedback bit is set depending upon which of 
the two estimates indicates a greater power. In other words, the receiver 80 
determines the power of the received pilot in the odd and even time slots, and 
determines whether the received power in the odd or the even slot was greater. 

FIGURE 6 is a block diagram of a preferred embodiment of a receiver (e.g., an MS) 
229 constructed according to this invention. FIGURE 6 illustrates the efficient use of 
the pilot format to both aid antennae adaptation and to demodulate incoming data. The 
pilot is used to extract a preliminary channel estimate c using the multiplier 202 and 
the accumulator 204. A sliding window filter 212 maintains the data channel estimate 
I. The preliminary data symbol estimates are generated by multiplying the data 
despreading sequence with the received signal using a multiplier 216 and 
accumulating the result in an accumulator 218. These preliminary symbol estimates 
are then delayed in a deskew buffer 219. The data channel estimate c is conjugated in 
a conjugation unit 214 and multiplied with the preliminary data symbol estimates by a 
multiplier 220 to generate the final data symbol estimates. A channel difference 
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estimate is simultaneously maintained through a separate data path. The preliminary 
channel estimate c is inverted for odd time slots and non-inverted for even time slots 
through the selection of the multiplexer 206 and application of ±1 in a multiplier 207, 
with this result accumulated in an accumulator 208 to generate c diff . When the 
measurement interval is complete, even and odd channel estimates are extracted from 
c diff and c using adders 210 and 224. The more powerful channel estimate is then 
selected using a decision block 228. 

In the case of resolvable multipath, the MS tracks several versions of the received 
pilot, making channel estimates for each one. In order to properly determine which 
Tx weight vector yields the best receive power, the MS additively combines the 
channel estimate powers from each path before performing a decision comparison. 
Similarly, if the MS 229 has multiple antennae, the decision is based on the sum of 
channel estimate powers over all antennae. 

The same method can be applied to soft handoff, where multiple base stations may be 
transmitting the same data to the MS with individually adjusted TxAA weight vectors. 
That is, the processing of the MS can be applied to each BS individually, and the 
odd/even power summations used to determine the feedback bit is the summation of 
the powers for each path and for each BS. All BS's receive the same feedback bit and 
adjust accordingly. This is almost effectively considering the weight vector to be a 
concatenation of the individual BS weight vectors. Each BS individually normalizes 
its weight vector. The described method provides a very simple manner for the TxAA 
to be distributed across several BS's with no increase in MS complexity. 

The following describes in greater detail the operation of the transmitter 50 and 
receiver 229 described above with reference to FIGURES 4 and 6, respectively. The 
following definitions are used in the description below with regard to the transmit 
adaptive antenna weighting technique: 

n Time index, at the nyquist rate (chip rate for DS-CDMA). 
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m Time index at decimated (slot) rate, m = \n I N\ . 

K Number of users to which the BS is transmitting. 

N Duration of the pilot even/odd banking slot. 

M Number of slots between weight updates. 

5 u(n) Transmitted vector at time n, [vl\ is the i th antenna. 

w(n) Antenna complex weights for the specific MS described, data channel. 

w even (n) Antenna complex weights for the specific MS described, even pilot 
channel. 

w odd (n) Antenna complex weights for the specific MS described, odd pilot 

10 channel 

x(n) Information bearing modulated signal for the specific MS described. 

p(n) Dedicated pilot for the specific MS described, 

gpiiot/traffk G a i n f° r pilot/traffic channel for the specific MS described. 

15 The complex baseband representation of the transmitted waveform at the antenna, 

"a 

Fy characterizing the amplitude and phase of the modulated radio carrier transmitted by 

the antenna (such complex baseband representations are well known, see, e.g., John 

ft;;; 

Proakis, Digital Communications, 1995, Third Edition, McGraw-Hill Inc., New York, 

fy 

I'D NY), is given as follows. The equation adds a subscript k differentiating the signals 

j. i 2 0 for the plural mobiles. 



j.4 



K-\l 



k=0\ 

+ ^(l«/^Jmod2 = 0)w Mv . *(») + (L«/^Jmod2 = l)w* 10 * *(n))p t (n) 

J 



For conciseness of presentation, the following text focuses on the present invention's 
25 method from the point of view of a specific mobile with an arbitrary index k, and the 
subscript k is omitted. 
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In particular, the odd transmit weight is applied to the dedicated pilot on odd time 
slots and vice versa for the even weights. To ensure proper coherent demodulation, 
the even/odd weights are preferably constrained by the following relationship: 



,v(;,)_ w ^^ ) + w ^ (n) 



It is worth noting that when the odd and even test weights are generated by the 
perturbation technique described below in more detail, then it is valid to consider that 



^ base\ n ) — 2 



and hence one may in that case use w base {ri) as the weight vector applied to the data 
channel. 

FIGURE 7 is a flow chart that outlines pertinent features of the operation of the 
transmitter 50. In block 100, a vector w base is initialized to an arbitrary vector with a 
unit norm. In block 102, a new test perturbation vector v(n) is generated. Each vector 
entry is a complex number and corresponds to one of the antennae 52 or 54; each such 
entry therefore corresponds to a test weight to be applied to the transmission of the 
pilot signal over a particular antenna. The test perturbation vector is preferably 
determined in a manner such that its entries appear mutually independent. As a point 
of reference, it is assumed that the expected value (or long term average over many 
realizations) of the magnitude squared of each entry of v is 2, so that 
£j|v| 2 )= 2 -(number of antennae). In this way, the amplitude of the applied 
perturbation is encompassed in the algorithm parameter f3 ', as described below. For 
example, v(n) could equal [±l±j, ... ±l±j] T (v(n) would have just two entries for the 
example transmitter of FIGURE 4), where the +*s and - s for the entries are randomly 
determined. However, the preferred v(n) is a complex Gaussian vector of 
uncorrected entries with variance of 2. 
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In block 102, the current value of v is stored as v old and a new value is determined for 
v. The old value is stored so that when the feedback is received at a later time (block 
1 16) the transmitter 50 knows which perturbation vector the feedback corresponded to 
and can adapt accordingly. In block 104, even and odd weights and data channel 
weights are updated based on the new value for v. This is summarized as follows. 
When beginning of test interval, 



v <= normalized test perturbation vector, new value every NM time indices; 
... w tee +^-v _ 



w 

f T even 



IK se +/MI 



w 



W even +yf odd 



2 
End. 

Note: the base weight vector is only updated when feedback is received. 

The normalization of the even/odd weight vectors is applied (the 2-norm is used) so 
that the total power transmitted in even and odd time slots is equal. This normalization 
prevents the MS 229 from simply selecting the weight vector with the larger transmit 
power. Rather, the MS 229 must select the weight vector which delivers the 
maximum power to the mobile for the given transmit power, which is determined by 
the traffic and pilot channel gains and g pilot . 

In block 106, the weights determined in block 104 are applied to the pilot signal. In 
particular, w even is applied to the pilot signal during even time slots and w odd is applied 
to the pilot signal during odd time slots, using the multiplexers 68 and 70 as is shown 
in FIGURE 4. In addition, w is applied to the traffic signal during both even and odd 
slots. The pilot signal and the traffic signal, as modified according to block 104, are 
transmitted from the antennae 52 and 54. More precisely, a sinusoidal carrier wave is 
generated for each of the antennae 52 and 54. Each of these carrier waves is 
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modulated by the traffic and pilot signals (modulating signals) as modified according 
to block 104. Viewed at from another perspective, the weight w adjusts the amplitude 
and phase of the plurality of radio carrier signals as those carrier signals are modulated 
by the traffic and pilot signals (modulating signals). 

In block 116, which is shown as disjoint from blocks 100-106 for reasons that will be 
described below, the transmitter 50 receives a feedback bit transmitted by the MS 229. 
This feedback bit is based on the decision from the MS 229, which indicates which of 
the even or odd channels resulted in a greater received signal power and determines 
the weight updates for the transmitter 50. In block 118, the transmitter determines 
whether the even weight yields greater power by checking the value of the feedback 
bit. In one embodiment, a value of "0" indicates that the even channel is better. If so, 
control passes to block 120, which sets equal to w even . In broader terms, w base is a 
function of w even and w odd . In the preferred embodiment, the function is an 'or' type 
function but more complex functions of weight vectors, including combinations 
thereof, are contemplated by the present invention. Otherwise, control passes to block 
122, which sets w base equal to w^. In block 124, the even and odd weights w even and 
and the data weight w are updated based on the new w^. This is summarized by 
the following pseudo-code: 
when feedback received 

if (feedback=0, indicating even channel is better) 
w toe +l-v 0 „ , 

base Ik^-vji' 

else 

w ■ ™ base -P 'Void . 

base K.-/?-v*r 

end 
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W = 



^ even ^ ^ odd . 



As noted above, block 116 (and related blocks) are shown as disjoint from blocks 100- 
5 106. This is because the processes represented by these two sets of blocks occur at 
different times, as can be seen from FIGURE 10, a timing diagram that shows a 
possible timeline for weight adjustment and channel measurement. Waveform 130 
shows the timing of updates to w base , which are performed in blocks 120 and 122. 
Waveform 132 shows updates to v, which are performed in block 102. Waveforms 
10 134 and 136 show updates to w even and w odd , which are performed in blocks 104 and 

124. Waveform 138 shows updates to the data channel weight vector w, which are 
performed in blocks 104 and 124. 

As shown by the waveform 132, the test vector v (step 102) is periodically generated 
15 independently of feedback received from the MS 229. The updates to w even and w odd 

caused by a new value of v are indicated by solid vertical lines in the FIGURE 10 
while the updates to w even and w odd caused by feedback received from the MS 229 are 
indicated by the dashed vertical lines in the Figure. 

2 0 The parameter /? is an algorithm constant which is to be selected to improve 
performance. A large p allows the weight vector to adapt more quickly, but also 
introduces additional noise on the weight vector applied as a larger adjustment is 
made. For example, when the weights are settled near their optimal value, each 
weight is still constantly adjusted by ±/?v % Also, ft determines the size of the test 

2 5 perturbation applied to the weight vector, and hence will effect both the ability of the 

MS to correctly determine which test weight vector provides the larger power, and the 
amount of interference seen by other MS's. A larger J3 will cause the difference 
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between the two test weight vectors to be larger, and hence it is more likely that this 
difference will exceed any noise contributions to the measurement by the receiver 80 
or receiver 229. A larger p also causes more "splatter", wherein the antenna pilot 
transmission test weights are perturbed more from the base weights, which are 
adjusted towards an optimal value, and hence the odd/even pilots cause more 
interference to other MS's than the data weight. As an enhancement to the present 
invention, it is possible to individually tune these two effects (i.e., the adaptation rate 
vs. weight noise tradeoff, and the noise immunity vs. splatter tradeoff) with a p x and 
a p 2 , as is described later. 

FIGURE 8 is a flow chart that outlines pertinent features of the operation of the 
receiver of FIGURE 5. In block 172 the even channel estimate is accumulated, which 
continues during the entire even time slot according to the decision determined at 
decision block 174. At the beginning of the odd time slot, the odd channel estimate is 
accumulated in block 176, which continues during the entire odd time slot according 
to the decision determined at a decision block 178. The even/odd channel 
accumulations continue until it is determined by block 180 that the measurement 
period is complete. When the measurement period is determined to be completed, the 
feedback decision is calculated and the resultant bit stored for transmission at the 
appropriate time in block 182. The odd and even channel estimates are then cleared in 
block 190 and the process is repeated for the new measurement period by returning to 
the processing block 172. At the appropriate time, which is independent of the timing 
of the channel estimation procedure just described, the MS starts transmission of the 
feedback bit, which is shown in processing blocks 184, 186, and 188. 

FIGURE 9 is a flow chart which outlines the pertinent features of the operation of the 
receiver of the preferred embodiment of FIGURE 6. In block 232 the sliding window 
data channel estimate is updated. The duration of the sliding window is shown as M 
slots, covering the duration of a measurement interval, as is preferred. The data 
channel estimate is applied to the preliminary symbol estimate in block 233. The state 
of the slot being processed is used in block 234 to determine the sign of the channel 
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difference accumulation by adding the preliminary channel estimate in block 248 if it 
is an even slot, or subtracting the preliminary channel estimate in block 250 if it is an 
odd slot. Block 240 determines if the measurement interval is complete, continuing 
the process by returning to block 232 if the interval is not complete. When the 
measurement period is determined to be completed by block 240, the feedback 
decision is calculated and the resultant bit stored for transmission at the appropriate 
time in block 238. The difference estimate is then cleared in block 236 and the 
process returns to block 232. At the appropriate time, which is independent of the 
timing of the channel estimation procedure just described, the MS 229 starts 
transmission of the feedback bit, which is shown in blocks 242, 244 and 246. 

FIGURE 10 shows a time line of this process with certain parameters selected by way 
of example. The accumulations for this example are over 12 time slots, the feedback 
bit is transmitted one time slot later, and the transmit weight update occurs as soon as 
the transmitter 50 receives it. In this case, the transmitter 50 updates the test 
perturbation vector v prior to receiving the update, so that the measurement of the new 
perturbation can begin. Because each update is fairly small, the measurement of the 
new perturbation will be insignificantly modified by the overlap. 

As previously described with reference to FIGURE 7, at the beginning of each test 
period, the transmitter 50 determines a new perturbation vector v and applies it to the 
old w base . As soon as the feedback arrives, this same new v is applied to the new w base . 
This allows for maximizing the measurement interval without delaying the next test 
perturbation. Because the updates are relatively small, the new perturbation will still 
yield a valid measurement result even though it is applied to both the old and the new 

W base- 

The invention may be implemented for "soft handoff ' systems. Soft handoff in DS- 

CDMA systems is an operational mode wherein multiple BS's, each using a different 

spreading code, transmit the same data to a single MS. The reception at the mobile of 

these multiple signals benefits performance by providing diversity. According to one 

embodiment of the present invention, soft handoff may be implemented in the 
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following manner. Each BS independently generates test weight vectors w odd and w even 
through the independent generation of test perturbation vectors v. Each BS transmits 
the pilot and data as previously described. The MS makes a single decision based on 
the summation of the channel estimate powers for all BS's, and transmits a single 
feedback bit as previously described. Each BS independently implements the weight 
adaptation as previously described, without regard to the fact that it is in soft handoff 
mode. In summary, the MS uses the summation of all BS powers for the decision, and 
each BS behaves exactly as though it would even if not in soft handoff mode. 

The application of the system described herein to the IS-95 standard and its derivative 
cdma2000 most likely would use N=64 chips (chip rate 1.2288MHz, giving 52.083us 
time multiplexing) and 1.25ms decision intervals. The length of the total even/odd 
accumulation is implementation specific in the mobile, but probably 8N, 12N, 16N or 
24N (1/2 of that for each of the even/odd accumulations). Another possible 
alternative would be 0.625ms decision intervals, in which case the total even/odd 
accumulation would probably not exceed 12N. These timings fit well into the existing 
specification as 64 chips is the original symbol duration from IS-95 and still a 
common time unit, and 1.25ms is a "power control group", the timing upon which 
over-the-air closed loop power control updates take place, and is 1/16 of the most 
common frame length, 20ms. It will be appreciated that the present invention is not 
limited to any of the specific quantities described above. 
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Some Alternate Embodiments 

It will be appreciated that it is possible to use more than 2 time slots and a 
corresponding number of test weight vectors, which may be generated from a 
corresponding number of test perturbation vectors. Also, more than 1 bit can be used 
as feedback, corresponding to a greater number of time slots. 

The generation of the test perturbation v can be implemented in many ways. Such 
generation should be pseudo-random so that over the long term the summation Sw" 
approaches an equi-diagonal matrix. The random elements can be generated as some 
probability distribution other than the binary one described {e.g., uniform distribution, 
Gaussian distribution). Rather than pseudo-random generation, the vector can be 
generated by cycling through some fixed sequence. For example, the sequence 

[10 0 ..0] T 
[j0 0..0] T 
[0 1 0 ..0] T 
[0j0..0] T 

can be selected in one embodiment, to scan through each weight independently. This 
satisfies the requirement on the long term summation of w H . 

The value of ft can itself be adaptive. For fast varying channels, the value could be 
made larger, and for slower channels it could be made smaller. This would allow for 
some optimization of J3 with respect to channel variation rate. 

The weight update can be further parameterized, allowing independent optimization of 
the adaptation rate vs. weight noise (/?,) and the measurability in noise vs. splatter 
(/? 2 ) tradeoffs, fi ', is then applied to the update of w^, while J3 2 is applied to the test 
weight vectors w odd and w even . This can be performed in accordance with the following 
pseudo-code: 
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when feedback received or new test interval 
if (new test interval) 

Void <= V 

v normalized test perturbation function 

end 

if (feedback received) 

if (feedback==0, indicating the even channel was better) 

, W base + A V old . 

else 



^ base 



^base Pl^old . 



end 

end 



even ^ 



+ ™odd . 

2 



W eve« + ^odd . 

w c= 



end 



The parameters used in the time line of FIGURE 10 are exemplary only and can be 
easily modified. 

The exact nature of the feedback channel is not specified. It can be implemented as a 
puncturing onto the data channel as the power control bit is implemented in IS95, or 
onto the pilot channel as in IS-2000, or in some other manner. 



The MS can generate the channel estimations in a different manner than described. 
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The weight normalization described above is a total power normalization. Other 
normalizations can be chosen. 

Alternate treatments of soft handoff can be used, though the previously described 
method is preferred. Rather than using only 1 feedback bit, a feedback bit may be 
included for each such base station (or each such BS that has adaptive Tx antennae 
enabled) so that the mobile can send commands to each, or alternatively can time- 
multiplex the feedback bits between the various BS's. The latter would result in a 
significant loss of performance to accommodate the slower adaptation times. 

The test weights w odd and w even are described as being applied in separate time slots. It 
is possible to use other orthogonal modulation techniques such that they may be 
extracted. For example, rather than an odd and even time slot, odd and even 
orthogonal codes (perhaps from a set of Walsh codes) could be used. This is not 
desirable for the preferred embodiment but may be desirable in other embodiments. 

While the transmission medium is assumed to be the propagation of radio waves from 
transmitting antennae to a receiving antenna (or antennae), the invention is obviously 
extendable to the propagation of waves in other media. For example, the medium 
could be acoustic waves in a fluid or gas, from transmitting transducers to a receiving 
transducer (or transducers). 
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Channel Autocorrelation Tracking Technique 

The channel autocorrelation tracking technique described herein improves weight 
vector tracking in a wireless communication system by taking advantage of 
knowledge of the correlation between elements of the forward channel gain vector. 
5 The present invention is particularly well suited for communication systems that 
utilize transmit adaptive antenna techniques. The present invention improves weight 
vector tracking by utilizing one of several techniques described hereinbelow. Several 
embodiments of the present inventive channel autocorrelation tracking technique are 
described hereinbelow. The present invention can utilize one embodiment or any 
10 combination of embodiments to improve weight vector tracking in wireless 
communication systems. 

In one embodiment of the autocorrelation tracking technique, the present invention 
extracts a biased gradient estimate by utilizing a correlated perturbation vector 
15 probing technique, where the perturbation vector is generated with a autocorrelation 
derived from an estimate of the channel vector autocorrelation matrix. The 
autocorrelation with which the perturbation vector is generated is selected such that 
the stronger subspaces of the channel autocorrelation are emphasized in the 
perturbation probing and update. This is desirable because these subspaces are more 

2 0 likely to be occupied by the channel vector at any given time than weaker subspaces, 

and thus the gradient estimate is biased to these emphasized subspaces. Specific 
embodiments use a perturbation autocorrelation given by an exponentiation of the 
estimate of the channel matrix added to an equi-diagonal matrix (e.g., scaled identity 
matrix). 

25 

The inventive channel autocorrelation tracking technique extracts a coarse gradient 
estimate by applying perturbation vectors to a weight vector; wherein the method 
utilizes a gradient feedback algorithm for adaptation of transmit antenna baseband 
weights using 1-bit feedback. In accordance with the present invention, the channel 

3 0 autocorrelation tracking technique determines a class of perturbation vectors for 

improving tracking when the channel gain vector contains correlated elements. A 
basic binary gradient feedback algorithm of the above-described TxAA Weighting 
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Embodiment that is adapted for use with the present invention is now described in 
some detail. 



The following definitions are used in the description below. 

receive unit The unit which is receiving the signals which are sent according to the 
transmit weight adaptation algorithm, generally referred to as the 
mobile station. This unit is also actively transmitting, as at a minimum 
the feedback information must be transmitted to the "transmit" unit, 
and probably there is other reverse traffic and signaling. 

transmit unit The unit which is transmitting the signals which are sent according to 
the transmit weight adaptation algorithm, generally referred to as the 
base station. This unit is also actively receiving, as at a minimum the 
feedback information must be received from the "receive" unit, and 
probably there is other reverse traffic and signaling. 
I I Matrix or vector norm by some defined norm, most often the 2-norm. 

I || 2 Matrix or vector 2-norm. For a vector x, the two norm is ||x|| 2 = Vx"x . 

For a matrix A, the two norm is ||a|| 2 =_Jnwx(x"A"Ax). 

J | Matrix "trace norm". For a K-by-K matrix A, the trace norm is the sum 

of the absolute values of the diagonals divided by the number of such 

1 K ~\ , ■ 11 

diagonals, ||A|| /r = — £ • The matrices of interest are typically non- 

negative definite, so that the diagonals are non-negative and 

l K ~ ] 

Wlfr Lnon-negativedef ~ K S^'* ' 



n 


Time index, at the nyquist rate. 


K 


Number of users to which the BS is transmitting. 


L 


Number of transmit antennas. 


N 


Duration of the pilot even/odd banking slot. 


M 


Number of slots between weight updates. 


u(n) 


L-by-l vector, transmitted vector at time n, [u]j is the j th antenna. 
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C((n) L-by-1 vector, forward channel gain vector, complex gain from each 

antenna of the Tx unit to the antenna of the Rx unit. 

c r (n) L-by-1 vector, reverse channel gain vector, complex gain from the 

antenna of the Rx unit to each antenna of the Tx unit. 

w(n) L-by-1 vector, antenna complex weights for the MS data channel. 

w even (n)L-by-l vector, antenna complex weights for the MS, even pilot channel. 

w odd (n) L-by-1 vector, antenna complex weights for the MS, odd pilot channel. 

v(n) L-by-1 vector, test perturbation to the complex weight. 

x(n) Information bearing modulated signal. 

p(n) Dedicated pilot. 

gpiiot/traffic Gain for the pilot/traffic channel for user. 

The basic binary gradient feedback algorithm of the TxAA Weighting Embodiment 
that is adapted for use with the present invention incorporates a simple 
implementation of the TxAA Weighting Embodiment. In accordance with the basic 
binary gradient feedback algorithm, transmitted waveforms are described for only one 
user (e.g., a receiver), while multiple waveforms are summed for multiple receivers. 
The complex baseband representation for K users of a transmitted waveform given by 
the complex conjugate of u(n) is given as follows, where the subscript k is added to 
denote the distinct users: 



K-l( 

+ g*,p,J(L"/^J mod2==0 K^ *(n) + in/N\mod2 = \)* kjlM *(nj)p k (n) 

Since the processing at the transmitting unit and the receiver unit is the same for all K 
users and is independent of the user K-l users, the subscript k will be excluded from 
the following detailed description. 

Thus, given the forward channel gain vector c fi the waveform received at the antenna 
of the receive unit is as follows: 
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r n {n) = u T {n)c f {n) 

= v? H (n)c f {n)g trajric 

Maximizing the received power for the user for a constrained weight vector norm 
results in the following ratio: 



This ratio, w, is the weight vector to which the algorithm converges. 

In accordance with the binary gradient feedback algorithm, odd transmit weights are 
applied to dedicated pilots during odd time slots and even transmit weights are applied 
to dedicated pilots during even time slots. In a receiver coherent demodulation 
embodiment, the algorithm constrains even/odd weights according to the following 
relation: 

) 2 

Receivers select the better of the even/odd to generate feedback. The algorithm 
generates channel estimates for even and odd weights. The summation of even and 
odd weights generates data channel estimates that are utilized for coherent 
demodulation. The difference of the squares of even and odd weights can be utilized 
to generate a feedback bit. In one exemplary embodiment, a logical zero value is 
arbitrarily defined as a condition when the even channel yields greater power and 
conversely a logical one value is defined as a condition when the odd channel yields 
greater power. In a resolvable multipath embodiment, receivers track several versions 
of a received pilot and determine channel estimates for each version. Prior to a 
decision comparison, receivers additively combine channel estimate powers from each 
path to determine which Tx weight vector yields the best receive power. Receivers 
transmit the same feedback bit(s) to all base stations when in soft handoff modes of 
operation. Base stations utilize this feedback to generate new transmit weights 
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according to the following pseudo-code: 



• when feedback received 



= normalized test perturbation function, new value every NM time indices 



if (feedback==0) 



w 



NM 



base 



NM 



-1 



else 



w,. 



NM 



-1 



-1 



end 



w 



ATM 



/c,eve« II s 



w 

end 



m7 



odd 



NM 



W eve« + ^ odd 



Test perturbation vectors can be generated by a plurality of methods, however they are 
normalized so that /? captures the tracking rate. Utilizing a larger /? yields faster 
tracking with larger tracking noise. The operation of the channel autocorrelation 
technique is now described. 
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The channel autocorrelation tracking technique described hereinbelow utilizes a 
perturbation vector generated as a random gaussian vector having a perturbation 
vector autocorrelation, Ry, The following equation represents the perturbation vector 
autocorrelation of the basic binary gradient feedback algorithm: 
R v =e(w h )=2I 

When the forward channel gain elements are independent, identically distributed (iid), 
and thus uncorrected, the forward channel autocorrelation estimate can be represented 
by the following equation: 



When the forward channel vector has this form and is uncorrected, the uncorrelated 
perturbation vector is appropriate. However, when the forward channel vector does 
contain correlation, a modified perturbation vector can be applied using an estimate of 
the forward channel gain vector autocorrelation matrix. 

The forward channel vector autocorrelation matrix estimate can be assisted by a 
parameterized approach, where the angle of arrival and angular spread of the multiple 
reflection paths between the transmitter antennas and the receiver are estimated. 
These parameters can be plugged into a parameterized autocorrelation matrix formula 
(several are available in the literature, see e.g. K. I. Pedersen et al "A Stochastic 
Model of the Temporal and Azimuthal Dispersion Seen at the Base Station in Outdoor 
Propogation Environments," IEEE Transactions Vehicular Technology, Vol. 49, No. 
2, March 2000) to provide an estimate, which could be used to assist in the overall 
matrix estimate. This parameterized approach could be used by itself or averaged 
with the other forward and reverse autocorrelation matrix estimates (e.g., above- 
described estimates). Note that the term "angle of arrival" is generally used to 
describe the angle of a particular path between antennas regardless of whether a radio 
wave is "departing" (z.e., being transmitted) or "arriving" (i.e., being received). 
Knowledge of the geographical position of a mobile station can also aid in the overall 



r = e(c r c / )= a 2 l ; for some scalar a . 
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matrix estimate. Any of several well-known methods of determining the geographical 
position of a mobile station can be used with the present invention, for example, a 
Global Positioning System (GPS). 



In general, it is unlikely that forward channel gain elements are uncorrected. The 
present invention advantageously forces the perturbation vector to be larger in the 
subspaces in which the channel vector c is more likely to occupy so that unnecessary 
weight noise is not added into subspaces that are unlikely or impossible for c to 
occupy. Accordingly, a reasonable perturbation vector autocorrelation can be 
represented by the following Equation 1 : 



Equation (1) 



R v =2 



1^ 

M 



The embodiment of Equation 1 uses a normalization of R f . In one embodiment, the 
present inventive channel autocorrelation tracking technique utilizes a norm referred 
to as a trace norm. Advantageously, the trace norm is relatively easy to calculate. 
One skilled in the communication art shall recognize that any method of calculating 
the norm can be utilized with the present invention without departing from the scope 
or spirit of the present invention. The trace norm of a matrix is defined as the sum of 
the absolute values of the diagonal elements of the matrix divided by the number of 
such diagonal elements. 

Equation 1 allows the present inventive technique to provide improved tracking when 
the channel vector is correlated and normal tracking (i.e., default perturbation of a 
TxAA weighting technique) when the channel vector is uncorrected. In accordance 
with the present channel autocorrelation tracking technique, a is an algorithm 
parameter, wherein 0<a<l. The parameter a is applied as in Equation 1, which 
ensures that all modes of the system can be exercised even if that mode is small in the 
autocorrelation matrix. In practical applications, the autocorrelation matrix is an 
estimate and the exact value is not known, and thus the algorithm parameter, a, 
ensures that all subspaces are exercised. For example, if the estimate is from the Tx 
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weight vector, the initial autocorrelation matrix has a nullspace, and a=l, then the 
nullspace will never be exercised by the weights, and thus never added to the adapted 
autocorrelation matrix estimate. Similarly, the actual channel autocorrelation matrix 
may change slowly over time so that subspaces which were previously null may 
become active. Setting a<\ ensures that all subspaces are exercised and the 
autocorrelation matrix estimate can continue to evolve without excluding the future 
possibility of active subspaces that are currently null. 

The operation of the channel autocorrelation tracking technique is summarized in the 
following steps: 

(a) determining channel autocorrelation matrix estimates (Rf, Ry); 

(b) determining a perturbation vector autocorrelation (R^) based on the 
autocorrelation matrix estimate; and 

(c) determining a weight vector (w) based on the perturbation vector 
autocorrelation by a TxAA algorithm and feedback from the 
receiver. 

Three embodiments of the channel autocorrelation tracking technique are described 
hereinbelow. The first channel autocorrelation technique of the present invention is 
now described in some detail. 

Channel Autocorrelation Tracking Technique - A First Embodiment 

The first embodiment of the channel autocorrelation tracking technique of the present 

invention extracts a coarse gradient estimate by utilizing a correlated perturbation 

vector probing technique, where the perturbation vector is generated with a 

autocorrelation derived from an estimate of the channel vector autocorrelation matrix. 

The autocorrelation with which the perturbation vector is generated is selected such 

that the stronger subspaces of the channel autocorrelation are emphasized in the 

perturbation probing and update. Specifically, the first embodiment utilizes a TxAA 

algorithm to generate perturbation vectors by applying perturbation vector 

autocorrelation matrix estimations using both forward and reverse channel estimators. 

The perturbation vectors are utilized by the TxAA algorithm to generate weight 

vectors. An exemplary TxAA algorithm is described above with reference to the 
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TxAA weighting embodiment. One skilled in the communication art shall recognize 
that any TxAA algorithm can be utilized with the present channel autocorrelation 
technique without departing from the scope or spirit of the present invention. The 
present channel autocorrelation tracking technique advantageously improves tracking 
when a channel gain vector contains correlated elements. The operation of the first 
embodiment of the channel autocorrelation tracking technique is summarized in the 
following steps: 

(a) determine a channel vector autocorrelation matrix estimate; 

i. calculate a forward channel autocorrelation matrix estimate (Rf); 

ii. calculate a reverse channel autocorrelation matrix estimate (Rf); 

(b) determine a perturbation vector autocorrelation (RJ based on the forward and 
reverse channel autocorrelation matrix estimates; and 

(c) determine a weight vector (w) based on the perturbation vector autocorrelation 
by a TxAA algorithm. 

The first embodiment of the channel autocorrelation tracking technique is a linear 
approach that utilizes the forward and reverse channels to estimate the autocorrelation 
matrix and apply the estimate to Equation 1. In accordance with STEP (a)(ii) of the 
first embodiment, the reverse channel estimate is determined by estimating the reverse 
channel gain vector through standard system procedures. For example, the reverse 
channel gain vector, c r (n), can be determined by multiplying the Rx waveform at the 

transmitting unit (e.g. a base station) by a local version of the reverse channel pilot 
sequence conjugated and filtering the result as is well known to anyone skilled in the 
pilot assisted CDMA communication art. In accordance with STEP (a)(i) of the first 
embodiment, the forward channel estimate is determined directly from the weight 
vector as shown in the following equation: 
c f (n) = yv(n). 

The first embodiment of the channel autocorrelation tracking technique utilizes the 
reverse channel as a component of the estimate of the forward channel autocorrelation 
because the autocorrelation of the reverse channel vector may be similar to the 
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autocorrelation matrix of the forward channel vector. This relationship is shown in 
the following Equation 2: 

Equation (2) R^R,. 

When Equation 2 is not valid, the reverse channel should not be used as an 
approximation for the forward channel. Thus, a„ which weights the reverse channel 
component of the overall estimation as will be described below, should be set equal to 
zero as described further below. 

In accordance with STEPS (a)(i) and (a)(ii) of the first embodiment, the forward and 
reverse channel autocorrelation matrix estimates are given by the following Equations 
(3) and (4), respectively: 

Equation (3) R f (n) = f> 7 (k)c f (n - % H {n-k). 

Equation (4) (n) = £ h r (k)c r (n - k)z r H (n - k) . 

where h r () and h£) are causal estimation filters with unit DC gain so that 

00 00 

2X(*)= 2>/W= !• 

k=0 Jt=0 

One simple filter adapted for use with the present invention is an exponential filter 

that can be represented by the following equation: 

0 k<0 
l(l-A)A k k>0 

Thus, yielding a simple recursive update that is represented by the following equation: 
R(n) = (l - zy(ny{n)+ JLR(n - 1) . 



h{k) = 



In accordance with STEP (b) of the first embodiment, the autocorrelation matrix of the 
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applied perturbation vector is given by the following equation: 



R v =2 



where a^and a r are algorithm parameters, wherein 0<a f <l, 0<a<l, and 0<(a f + a r )<l. 

The algorithm parameters, a f md a r , can be determined in any convenient manner such 
as an empirical method {e.g., based on field measurement and calibration of antenna 
units). In the first embodiment, the approximation referenced in Equation 2 can be 
utilized when the basestation uses the same antennas for Tx and Rx. Also, depending 
on the geometry and radio carrier frequencies used, the reverse channel may or may 
not be useful in generating the perturbation autocorrelation matrix. If the reverse 
channel is not useful then a, is set to equal zero. 

In accordance with the first embodiment, the time constants of the autocorrelation 
matrix estimations would be much larger than the time constant of the tracking given 
by In one exemplary embodiment, the autocorrelation matrix estimations are 
many times longer in duration (e.g., 100 to 100,000 times longer). 

The STEP (c) of the first embodiment can be performed by any TxAA algorithm. An 
exemplary TxAA algorithm that can be utilized with the first embodiment of the 
present invention is described above with reference to the TxAA weighting 
embodiment. 



The first embodiment of the channel autocorrelation tracking technique is a linear 
approach that utilizes the forward and reverse channels to estimate the autocorrelation 
matrix and apply the estimate to Equation 1. One skilled in the communication art 
shall recognize that other approaches such as non-linear approaches can be utilized 
with the present invention without departing from the scope or spirit of the present 
invention. One such non-linear approach is now described with reference to the 
second embodiment of the channel autocorrelation tracking technique. 
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Channel Autocorrelation Tracking Technique - A Second Embodiment 
The second embodiment of the channel autocorrelation tracking technique of the 
present invention extracts a coarse gradient estimate in a similar manner to the first 
embodiment. Specifically, the second embodiment utilizes eigendecompositions to 
generate non-linear functions of the autocorrelation matrix. These non-linear 
functions can provide improved tracking by 11 expanding" and "compressing" the 
tracking rate within the different subspaces of the autocorrelation matrix estimate. 
Expanding the tracking rate is defined as emphasizing stronger eigenmodes and de- 
emphasizes weaker eigenmodes. Compressing the tracking rate is defined as bringing 
the tracking rates of strong and weak eigenmodes closer. 

As described further below, a parameter p 9 applied as a matrix exponent, determines 
the degree of expansion or compression of the tracking rate. In the limit when p=0, 
the autocorrelation matrix is not used to assist tracking and all subspaces are tracked at 
the same rate. In the limit when /?-»ao, only the dominant subspace is tracked and the 
others are ignored. When p=\ the second embodiment operates in a similar manner to 
the first embodiment. The operation of the second embodiment of the channel 
autocorrelation tracking technique is summarized in the following steps: 

(a) determine a perturbation vector autocorrelation matrix estimate (Rf, based 

A p A p 

on non-linear channel subspace weighting matrices {e.g., R/ ,Rr and 
& f f ) utilizing eigendecompositions; 

(b) determine a perturbation vector autocorrelation (RJ based on the forward and 
reverse channel subspace weighting matrix; and 

(c) determine a weight vector (w) based on the perturbation vector autocorrelation 
by a TxAA algorithm. 

The second embodiment of the channel autocorrelation tracking technique is a non- 
linear approach that can improve performance as the technique expands or compresses 
the tracking rate versus noise tradeoff for the various eigenmodes of the 
autocorrelation matrix. The second embodiment utilizes a tracking rate ratio, y , to 
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expand or compress the tracking rates. The tracking rate ratio is defined as the ratio of 
one eigenmode to another eigenmode. Two exemplary methods of the second 
embodiment of the channel autocorrelation tracking technique are described 
hereinbelow. 

In accordance with STEPS (a)(i) and (a)(ii) of the second embodiment of the channel 
autocorrelation tracking technique, Equations (3) and (4) described above with regard 
to the first embodiment undergo eigendecompositions. These eigendecompositions 
can be represented by the following equations: 

R f =Q f A f Q»; 

R,=QAQf; 

where the matrices Q are comprised of the eigenvectors of the 
corresponding matrices R; and 

the matrices A are diagonal, containing the eigenvalues of R. 

These eigendecompositions can be performed by any standard method. One such 
method is described in a book by G.H. Golub and C.F. van Loan, entitled "Matrix 
Computations, Third Edition", published by Johns Hopkins University Press in 1996, 
which is hereby incorporated by reference for its teachings on eigendecompositions 
and matrix algebra. 

In accordance with STEP (a) of the first exemplary method of the second 
embodiment, these eigendecompositions are used to generate non-linear functions of 
the autocorrelation matrix that can be represented by the following equations: 

R/ =Q r A/Qf 

When p is an integer the power operation can be performed without the 

eigendecomposition, however, performing the eigendecomposition is effective for 

generating the correlated perturbation vector because the method can simply use the 

eigendecomposition as the orthogonal basis of the random vector generation. The 
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selection of algorithm for performing the eigendecomposition is dependent on the 
value of p and the number of dimensions to the matrices. 

In accordance with STEP (b) of the first exemplary method of the second 
embodiment, the autocorrelation matrix of the generated perturbation vector can be 
represented by the following Equation 5: 



Equation (5) R v =2 



, M l R ' 

where a f and a r are algorithm parameters, wherein 0<a f , a,<l and 
0<(a f + a r )<L 

The STEP (c) of the second embodiment can be performed by any TxAA algorithm. 
An exemplary TxAA algorithm that can be utilized with the second embodiment of 
the present invention is described above with reference to the TxAA weighting 
embodiment. An exemplary second embodiment of the channel autocorrelation 
tracking technique has been described. One skilled in the communication art shall 
recognize that alternative methods of performing the second embodiment can be 
utilized with the present invention without departing from the scope or spirit of the 
present invention. A second exemplary method of the second embodiment is now 
described. 



The second exemplary method of the second embodiment of the present invention 
performs the operation STEPS (a) and (b) in accordance with the following 
description. STEP (c) is performed in a similar manner to the first exemplary method 
of the second embodiment. 



In STEP (a) of the second exemplary method of the second embodiment, the method 
utilizes a joint exponentiation of forward and reverse estimate components, according 
to the following equation: 
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R 



R, 



/ R, 
R. 



= QfAf/Qf„- 



In STEP (b) of the second exemplary method of the second embodiment, the method 
utilizes the following Equation 6: 



Equation (6) 



R 

R.,=2 {a f + a) ¥ ^ 7 r^-a r -a f 
R 1 



The present invention can use either exemplary methods of the second embodiment to 
improve transmit antenna weight tracking in a communication system. The present 
invention may use the second exemplary method instead of the first exemplary 
method based on computational complexity. For example, Equation 6 has less 
computational complexity than Equation 5. 

The third exemplary method of the second embodiment of the present invention 
performs the operation STEP (b) in accordance with the following description. 
STEPS (a) and (c) are performed in a similar manner to the first exemplary method of 
the second embodiment. 



In STEP (b) of the third exemplary method of the second embodiment, the method 

utilizes the following equation: 

f * Y 
a, 



R v =2 





• + a r 1. II 


||R f 


R r 



The third exemplary method is advantageously less complex than the first exemplary 
method. To wit, only one eigendecomposition must be calculated versus two. A 
further advantage is the inclusion of the identity matrix prior to the 
eigendecomposition, which ensures a numerically robust decomposition. Such good 
conditioning for robust decomposition may not be available from the raw values of R 
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without the additional diagonal weighting, because the matrices R may not be well 
conditioned. An exemplary implementation of the channel autocorrelation tracking 
technique is now described. 

An Exemplary Implementation of the Channel Autocorrelation Tracking Technique 
The exemplary implementation of the channel autocorrelation tracking technique 
identifies exemplary values for parameters such as ft and p. These parameters can be 
determined from simulations and field data. In simulations of bit error rate for 
tracking of uncorrected fading with systems having L=3 transmit antennas, the 
exemplary implementation sets /? equal to a value of O.L The exemplary 
implementation sets p equal to a value of 2 because it is desirable to leave the "all- 
mode" excitation at 1/10 of the normal excitation as this induces very small excess 
weight noise, and thus it is advantageous to set p greater than 1 to take advantage of 
highly correlated environments. Thus, the exemplary implementation comprises the 
following algorithm parameters for L=3 antennas: 



The exemplary implementation utilizes the non-linear implementation of the second 
embodiment of the channel autocorrelation tracking technique (i.e., Equation 6). If 
Equation 2 is valid (i.e., the reverse channel autocorrelation is useful in determining 
the forward channel autocorrelation), then a f =a=0A5, otherwise a r =0 and ^=0.9, In 
general, Equation 2 will be invalid when the antennas are spaced far apart or the 
duplex frequency is large. 

To facilitate a better understanding of the usefulness of the reverse channel 
autocorrelation in determining the estimation, an example is now described. A system 
having an algorithm measurement feedback period of 0.625ms requires 10,000 sample 
subspace weighting matrices with a/=0, and thus requires a autocorrelation estimation 
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filter time constant of 6.25 seconds. If the reverse channel can be used to help 
determine the estimation, then 10,000 samples (half from the forward channel and half 
from reverse channel) requires a autocorrelation filter time constant of 3.125 seconds 
on each of the forward and reverse estimators. Thus, the reverse channel 
autocorrelation is advantageously utilized to determine the estimation. 

FIGURE 1 1 is a block diagram of an exemplary apparatus adapted for use with the 
first embodiment of the channel autocorrelation tracking technique. FIGURE 1 1 is a 
baseband equivalent model of the first embodiment. The exemplary apparatus 
generates perturbation vectors without subspace emphasis through matrix 
exponentiation by the exponent p. The following description regarding FIGURE 1 1 
describes the operation for a single transmission x(ri), however, the present invention 
is not limited to a single transmission. Those skilled in the art shall recognize that the 
present technique can operate in a parallel manner using multiple transmissions to 
multiple receivers without departing from the scope or spirit of the present invention. 

As shown in FIGURE 11, the exemplary transmitter unit apparatus comprises an 
antenna array 508, a vector duplexer 510, an estimation filter 524, a reverse link 
subspace weighting matrix generator 520, a forward link subspace weighting matrix 
generator 518, a perturbation vector generator 530, multipliers 504, 512, 522, 536, 
538, normalizes 526, 528 and adders 506, 534, 532, The antenna array 508 
comprises multiple antennas, where each antenna is used for both transmission and 
reception. The vector duplexer 510 represents a single duplexer for each antenna of 
the antenna array 508 so that transmission and reception vectors are isolated. The 
antenna array 508 inputs a received vector to the vector duplexer 510. The vector 
duplexer 510 inputs the received vector to the multiplier 522. 

The multiplier 522 multiplies the received vector by a scalar reverse link pilot 
sequence 552. The multiplier 522 inputs a resultant vector to the estimation filter 524. 
The estimation filter 524 generates a received channel vector estimate c r (rc) and 
inputs the received channel vector estimate c r (n) to the reverse link subspace 
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weighting matrix generator 520. The reverse link subspace weighting matrix 
generator 520 generates a reverse subspace weighting matrix. 

A forward weight vector 514 inputs to the multiplier 504 and is multiplied by a scalar 
traffic sequence 500 to produce a forward resultant vector. The forward link subspace 
weighting matrix generator 518 inputs the forward weight vector 514 to generate a 
forward subspace weighting matrix. The reverse subspace weighting matrix is 
normalized and scaled by the normalizer 528 and multiplier 538, respectively. The 
forward subspace weighting matrix is normalized and scaled by the normalizer 526 
and multiplier 536, respectively. The outputs of the multipliers 538, 536 are added to 
a scaled identity matrix through adders 534 and 532. The adder 532 inputs a resultant 
matrix 1^ to the perturbation vector generator 530 whereat the resultant matrix Ry is 
used as an autocorrelation matrix of a generated vector v. This generated vector v is 
used to generate even/odd weight vectors 516 which are applied through multiplier 
512 to the forward pilot sequence 502. The generated vector v is also used in updating 
the traffic weight vector w 5 14 with feedback 540 from the receiver. 

Summary 

The present channel autocorrelation tracking technique improves weight vector 
tracking in a wireless communication system by taking advantage of the observation 
that tracking can be improved when the channel vector representing the gain from 
multiple transmit antennas to a single receive antenna is comprised of correlated 
elements. In a first embodiment of the autocorrelation tracking technique, the present 
invention extracts a coarse gradient estimate by utilizing a correlated perturbation 
vector, where the perturbation vector autocorrelation is determined from an estimate 
of the forward channel vector autocorrelation, and the extracted gradient is used to 
update the transmission weight vector. 

In a second embodiment of the channel autocorrelation tracking technique, the present 

invention adjusts the tracking rate of the weight vector in the different subspaces of 

the channel estimate autocorrelation matrix by emphasizing or de-emphasizing the 

differences of the eigenmodes of the channel vector autocorrelation. Specifically, the 
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autocorrelation applied to the perturbation vectors is an exponentiated version of the 
estimate of the channel autocorrelation. In one embodiment this exponentiation is 
performed using eigendecompositions. These non-linear functions can provide 
improved tracking by "expanding" and "compressing" the tracking rate of the different 
subspaces of the channel subspace weighting matrix. 

A number of embodiments of the present invention have been described. Neverthe- 
less, it will be understood that various modifications may be made without departing 
from the spirit and scope of the invention. For example, the weight vector update can 
be performed by the TxAA algorithm or another feedback algorithm that comprises 
one-bit feedback. 

Accordingly, it is to be understood that the invention is not to be limited by the 
specific illustrated embodiment, but only by the scope of the appended claims. 



